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Abstract 

O . We study high mass dilepton production from giuon-induced processes, gq — > 

b:;;! ', q^*, gq — > qj*, and gg — > qq'j*, in a thermally equilibrated but chemically 

f_^ I under-saturated partonic matter that is expected to be created in the initial 

Qh' stage of ultra-relativistic heavy ion collisions. Regulating the divergence in these 

■^ . processes by the thermal quark mass, we find that gluon-induced processes are 

ff^ I more important than the leading-order qq annihilation process as a result of the 

larger number of gluons than quarks in the partonic matter. The dependence of 
the thermal dilepton yield from the partonic stage of heavy ion collisions on the 



(N 
> 



C^ I initial conditions for the partonic matter is also studied. We further discuss the 

CS ' feasibility of observing thermal dileptons from the quark-gluon plasma in heavy 

f— ^ . ion experiments. 

^ '. 1 Introduction 






One of the signatures for the formation of a quark-gluon plasma (QGP) in heavy ion 
collisions at the Relativistic Heavy Ion Collider (RHIC) is the thermal dileptons it 
pj , produces. Many studies have been carried out to estimate the dilepton yield from 

the QGP. In earlier works, contributions from the leading-order (LO) qq annihilation 
process [0] and processes up to 0{as) p, |^ have been evaluated for a thermally and 



X 



;h ' chemically equilibrated QGP. Although thermal equilibrium can be reached in the 



c^ 



partonic matter created in heavy ion collisions at RHIC energies, chemical equilibrium 
is, however, unlikely to be established [^ |^. Later studies have thus included the 
effects of chemical non-equilibrium on thermal dilepton and photon production from 
QGPi,a|,i. 

Since the partonic matter is dominated by gluons, the thermal dilepton yield from 
the LO qq annihilation process is thus much smaller than that from the initial Drell- 
Yan process and final charm decays |T^. However, the dominance of gluons in the 



initial partonic matter provides the interesting possibility that dilepton production 
from gluon-induced processes may become more important than the LO qq annihilation 
process, even though they are suppressed by powers of as- In this paper, we shall study 
high mass thermal dilepton production from gluon-induced processes, i.e., the Compton 
processes involving one gluon in the initial state {gq -^ g7* and gq -^ q^y*) and the 
two-gluon fusion process {gg — > qq^y*)- In particular, we shall consider these processes 
only at leading-order in a^, i.e., at 0{as) and 0{a'j.) for the Compton and two-gluon 
fusion processes, respectively. 



This paper is organized as follows. In Section ^ we introduce the thermal quark 
mass in a non-equilibrium partonic matter at finite temperature, which will be used 
in Section |^ to regularize the divergent dilepton production cross sections for gluon- 
induced reactions. Thermal dilepton rates are then calculated in Section ^ for an 
under-saturated partonic matter in thermal equilibrium. In Section ^, we estimate 
the total thermal dilepton yield in heavy ion collisions at RHIC energies and study 
its dependence on the initial conditions for the partonic matter. We then compare it 
with the Drell-Yan yield to see if thermal dileptons from the QGP can be observed in 
experiments. Finally, a summary and discussions are given in Section ^. 

2 The Thermal Quark Mass 

For a partonic matter in thermal equilibrium at temperature T but out of chemical 
equilibrium, the parton distribution can be written in the Jiittner form, 

f. = h (I) 

where i denotes different parton species, /3 = 1/T, and u is the four velocity of the 
local comoving frame. In the limit of small fugacities (Aj ^ 1), the parton distribution 
functions reduce to the Boltzmann form, 

/, ~ X,e-^^-P. (2) 

Assuming that the partonic matter is both baryon and flavor symmetric with Nf quark 
flavors, i.e., Xq^ = Xq. = Xq {i = l,Nf), then two fugacity parameters, Xg and Xq, are 
needed to specify the degree of chemical equilibrium. 

Since chemical equilibration of the partonic matter created in heavy ion collisions 
at RHIC is likely to be slow, we shall use the Boltzmann distribution of Eq.(^ in the 
following calculations unless specifled otherwise. In this case, the parton density in the 
local comoving frame is directly related to the fugacity, 

. dnT^ , {da + dq)NfT^ 
% = \^^ ^<i + ^5 = \— 2 ' (3) 

if partons are taken to be massless. In the above, di is the degeneracy factor of parton 
2, i.e., dg = 2{N^ - 1) and dg = dg = 2N^. 

The thermal quark mass, denoted as /i, is given by 



Nl-l r d^p 1 



/x^ = Anas^-— —^,-ifg + fq), (4) 



where as = g'^/{4:7c) is the QCD strong coupling. For a QGP in chemical equilibrium, 
this gives a thermal quark mass /x^ = g'^T'^/6 when using the exact parton distribution 
of Eq.(|I]). For an under-saturated partonic matter described approximately by Eq.(^, 
the thermal quark mass becomes 
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l^' = ^!-^iXg + Xq). (5) 



3 Dilepton Production Cross Sections 

Leading-order processes in the strong coupling constant as for dilepton production from 
parton-parton interactions are qq annihilation {qq —>■ 7*), Compton processes [gq -^ 
57* and gq — *• 57*), two-gluon fusion {gg — > qq'j*), and two-quark {qq -^ qq'J*) and two- 
antiquark bremsstrahlung {qq — > qq'y*). They are of order 0{1), 0{as), 0{a1), 0{a1) 
and 0{a'^), respectively. In the present study, we neglect two-quark and two-antiquark 
bremsstrahlung as they do not involve gluons in the initial state and are also suppressed 
by a^ relative to the LO qq annihilation process. Furthermore, we neglect 3—^2 
processes as they are expected to be less important when parton fugacities are much 
below one. The processes included in our study for dilepton production are shown in 

Fig-E- 

The cross section for dilepton production from partons can be generally written as 






(6) 



In the above, indices i and j specify initial-state partons, s = {pi + pjY , and a is the 
fine structure constant. The quantity Wij is related to the square of the transition 
amplitude. For the leading-order qq annihilation process, it is simply 



W, 



m 



e^5(l-0, 



(7) 



with i = W^/s. 

For the Compton processes, the squared matrix element, after averaging over initial 
and summing over final spins and colors, is given by 



\M .|2 



\M - - *|2 



{AnYaasel 



t s 2M^{s + t-M^ 

-- - ^ + ^ . 

St st 



The total cross section, given by a = /(it|Mp/(167rs^), is divergent. Following 
previous studies of low mass dilepton and photon production from the QGP |TT, |T^, T3|, 



we use the thermal quark mass as the screening mass for regularizing all divergence in 
the cross sections. We thus replace 1/t with l/(t — /i^), and Eq.(||) then becomes 



I A/f 1 2 
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2M^{s + i- M^ 



(9) 



After integrating over t, we obtain 



W„„ = W, 
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where 



co(0 



1 + 2^ - 3^2 
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l-2{ l-{ + A_5 
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There are no studies on dilepton production from the two-gluon fusion process, 
99 ~^ QQl* 1 ^^ finite temperature. However, processes involving heavy quarks in the 
final state such as gg — > QQ'y* have been studied in determining the K-factor in Drell- 
Yan processes in the limit of M^ ^ ttiq [n, 15]. These results can be adopted for 
the production of dileptons with invariant mass M > 2 GeV. Since M^ ^ /i^ unless 
the temperature is extremely high, we can replace the heavy quark mass rriQ in Ref. 



15| with the thermal quark mass /i and obtain 



W, 
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Ci(01n2^ + C2(01n^ + C3(0 
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In the above, the functions ci(^), C2(0 and C3{^) can be obtained from Ref.|Tj], where 
the calculation was carried out in the MS scheme. 



4 Thermal Dilepton Rates 

For a partonic matter in thermal equilibrium at temperature T but out of chemical 
equilibrium, the thermal dilepton rate is given by 



dRij 
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dM 
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with Wrei = (pi ■ P2)/(yPiP2), Sgg = 1/2 accouutiug for the two identical gluons in the 
initial state, and Sij = 1 otherwise. The summation is over the parton spin and color 
indices. 

The integrals in the above equation can be evaluated using the rapidity and trans- 
verse momentum variables. 



s = M^, P^ = (M^o cosh Fo, Pi.0 cos $o, P±o sin $o, M^o sinh Yo) , 



Pi = p±i (cosh 2/1, cos 01, sin 01, sinhyi) 



(14) 



where M^o = ^M^ + Pl^. Inserting /c?^Po'^^(Po - Pi - P2) into Eq.(|g), the deha 
function can then be used to remove the / d^p2 and / dp±i integrals. After integrating 
over $0 and redefining yi — Yq and 0i — $0 as yi and (pi, we obtain 



dRij 
~dM 



dd.kXjS,, , ^MndY,dM 
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(Mj_o cosh yi - P^o cos c/^i^ 
Carrying out the integrals over Yq and 0i and then those over yi and M±o gives 
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Summing over quark flavors, we finally have 
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for A^ 
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For the strong coupling constant, we use the LO value with three quark flavors. 
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with A = 200 MeV and the normalization scale chosen to be the dilepton mass M. 

In the limit M^ S> /x^, the thermal dilepton rates depend on the coupling constant 
and fugacities as 

dRn 



'^qq 2 \ 2 
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In Fig. 1^, we show the thermal dilepton rates from a quark-gluon plasma in full 
thermal and chemical equilibrium at T = 570 MeV. It is seen that for dileptons with 
mass above 2 GeV the LO qq annihilation process gives the dominant contribution, 
while the Compton processes give a small contribution, and that from two-gluon fusion 
process is even smaller. These results are not unexpected as the Compton and two- 
gluon fusion processes are suppressed by as and a], respectively. 

The situation is, however, very different if partons are not in chemical equilibrium. 
As an example, we show in Fig. ^ the thermal dilepton rates for T = 570 MeV, 
Xg = 0.060 and Xq = 0.0072. The reason for choosing these values will be given in 
the next Section. In this case, the Compton processes produce more thermal dileptons 
than the LO qq annihilation, and even two-gluon fusion processes give a comparable 
contribution to lower- mass dileptons. As a result, the total thermal dilepton yield is 
significantly enhanced by gluonic processes. 

The appreciable difference between the results shown in Fig. |] and Fig. ^ can 
be qualitatively understood from the scaling behavior given in Eq.(^). Because of 



the quadratic dependence on parton fugacities, gluonic processes become dominant 
over the qq process when gluons greatly outnumber quarks {Xg/Xg ^ 1). Even if 
Xg/Xq is fixed but Ag ^ in an under-saturated partonic matter, both Rgg+gg/Rgg and 
Rgg/Rgg becouic logarithuiic divergent, indicating that gluonic processes become more 
important than qq annihilation when the parton system is further away from chemical 
equilibrium. 

To illustrate this quantitatively, we show in Fig. ^ the full-equilibrium results in 
Fig. multiplied by the quadratic fugacity factors. In this case, the Compton processes 
are already as important as the qq annihilation process. Further increase of the gluonic 
contribution when going from Fig. ^ to Fig. |^ is due to the logarithmic terms in 
Eq . (^2]) , which come from the thermal quark mass in the regularized partonic cross 
sections given by Eqs.(|TU|-|r^). 

We also see in Figs. |^-^ that the LO qq annihilation process is most effective in 
producing high invariant-mass dileptons because all incoming energy can be used. In 
contrast, there are one or two partons in the final state in the Compton and two- 
gluon fusion processes to carry away part of the incoming energy. To produce high 
mass dileptons from these processes thus requires more energetic partons, which are 
less available. Therefore, the three curves in Fig. |^ corresponding to the above three 
processes exhibit increasingly larger slopes. 

5 Thermal Dilepton Yield at RHIC 

In this Section, we estimate the thermal dilepton yield in central Au+Au collisions at 
RHIC by integrating the thermal dilepton rate over the space-time four volume, which 
will be taken from a schematic model. Also, we will study the feasibility of observing the 
thermal dilepton signal in experiments by comparing it with the Drell-Yan background. 
Starting from Eq.(0) and using d^x = rdrdrid'^xx, we first integrate over rj instead 
of Yq and then over (pi, yi and M±o as in previous section. This gives 

^ ^ /dr„d..dV*4^//MoM'rA,(M„/r)^(M,M.), (23) 

which then leads to 

^^^^- - f^' dTTd'xA{T). (24) 



dW^dYo Jro dM^ 

In the above, Tq is the initial proper time when the parton system first reaches thermal 
equilibrium, and Tf is the final proper time when the temperature drops to Tc (chosen 
to be 200 MeV) and the partonic matter makes a transition to the hadronic matter. 
The rate dRij/dM"^ depends on r implicitly through the time evolution of the variables 
Aj,T and /z. 

We note that while the pair rapidity Yq of the two partons in the initial state is 
equal to the pair rapidity Y of the dilepton in the LO qq annihilation process, they are 
different in gluon-induced processes. However, if we assume that the expansion of the 
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partonic matter follows the Bjorken scaling, the dependence of the thermal dilepton 
yield on Y and Yq must have the following functional form: 

dN 



F(r-yo,M), (25) 



dM^dYdYo 

as a result of boost invariance of the integrand on the right-hand side of Eq. (|T3|) . This 
gives 

dN dN 



dM^dY dM^dYo 
The same argument applies to r], and it is easy to verify that 

dN dN 



(26) 



(27) 



dM^dr] dM^dYo' 

5.1 Time evolution of parton fugacities and temperature 

First, we need to specify the time evolution of partons in the quark-gluon phase created 
in heavy ion collisions at RHIC Earlier studies |]lOl based on kinetic equations including 



2 <-^ 3 gluonic processes have shown that thermal equilibrium is achieved first, and both 
gluon and quark fugacities increase as the system evolves. However, since the quark- 
gluon plasma has only a finite lifetime of a few fm/c, chemical equilibrium is unlikely 
to be established, and parton fugacities thus remain well below the value of one. To 
take into account such effects, we follow Fig. 25 of Ref.[|I^ by simply assuming that 
the parton fugacities increase linearly with the proper time, 

A, = Xlo (:^) • (28) 

In an ideal hydrodynamic description, temperature of the partonic matter decreases 
as T oc 1/r^/^, which then leads to a decrease of the rapidity density of transverse 
energy as dE^/dy oc tT^ oc 1/r^/^. In the parton cascade model, the time evolution 
of the rapidity density of transverse energy depends sensitively on the magnitude of 



parton cross sections and densities |T^ . Since we do not expect the rapidity density of 



transverse energy to change more than 25% during the partonic stage of RHIC heavy 



ion collisions |T^, |T^, we assume for simplicity that dEn^/dy, given by 



^ = ^<P^> oc (e ^^.o) tT'T oc {rT^y , (29) 

is a constant. This then leads to the time evolution of temperature as 

T = To (^^ j . (30) 

In the above, we have used 

— = TonR^2^ni (31) 



at central rapidity at tq and Ra — 1.12 A^^^ fm. Also, we neglect the transverse 
expansion throughout this study. The effect of transverse expansion on the production 
of dileptons of mass above 2 GeV has been shown to be very small in heavy ion 



collisions at RHIC energies |]T8|. This is expected because high mass dileptons are 
mainly produced in the early stage of the partonic matter, when the transverse flow 
has not much developed. 

In Fig. ^ the time evolution of fugacities and temperature given in Eqs. (|2^) and 
(pop are shown. We see that the temperature decreases faster than that given by the 
scenario of entropy conservation, i.e., T oc 1/r^/^, as parton production increases the 
entropy. 

5.2 The default scenario 



As in Ref. |T0[, where the parton equilibration is studied using kinetic equations, we 
assume that at the proper time ro=0.7 fm/c, thermal equilibrium is established at a 
temperature Tq = 570 MeV. We also take the same gluon and quark densities at tq as 



in Ref. [110], and obtain the following values for the initial parton fugacities |]T^ 



Ag,o = 0.060 and Ag,o = 0.0072. (32) 

With these values, we find that quarks and antiquarks make about 15% of the partons 
and the thermal quark mass is about 75 MeV if the strong coupling constant is evaluated 
at the normalization scale of 2 GeV. Subsequent expansion of the partonic matter is 
assumed to follow the Bjorken scaling. 

We have evaluated the thermal dilepton yield using Eq.(^), and the results are 
shown in Fig. |^. Also shown is the yield from primary Drell-Yan hard collisions using 
a K-factor of 1.7, obtained from fitting the E772 data, to take into account higher- 
order contributions. As in thermal dilepton rates shown in Fig.^ we find that the 
Compton processes produce more dileptons than both the LO qq annihilation and the 
two-gluon fusion processes. This also indicates that the conventional way of using a K 
factor of 2 to account for all higher-order contributions |§] underestimates the possible 
large enhancement from gluonic processes. We note that, despite the enhancement 
from gluonic processes, the total thermal dilepton yield is still a factor of 30 or more 
(depending on the dilepton mass) below the Drell-Yan yield. It is thus unlikely that 
thermal dileptons from the partonic matter can be observed in heavy ion collisions at 
RHIC. 

5.3 An upper limit 

The small thermal dilepton yield discussed in the above default scenario is mainly due 
to the small number of quarks and antiquarks compared to that of gluons assumed in 
the initial conditions for the partonic matter. While we have used the ratio of initial 
quarks to gluons that is consistent with minijet calculations based on perturbative 
QCD together with initial and final state radiations, it is possible that additional soft 
partons may be produced. Also, there may be already some chemical equilibration 
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before the thermal stage starts at tq. Both could then lead to initial quark fugacity 
that is larger than the estimate from minijet calculations. 

To obtain an upper limit for the thermal dilepton yield, we assume that quarks and 
antiquarks are in relative chemical equilibrium with gluons at the initial time tq, i.e., 
Xq = Xg, although both are still out of full chemical equilibrium. Requiring the same 
value of total dET/dy at central rapidity as before and using the same tq and To, we 
obtain from Eq. (p9[) the following initial parton fugacities for the upper limit (UL): 



it/L _ \UL _ dg Xgfi + {dq + dq)Nf Xqfl 
dg + {dq + dq)Nf 



\UL _ ^UL _ "^g^gfl "T l"-g "T u,qji^f/\qfi 



This gives a ratio of 3/5 for the sum of quarks and antiquarks to the total number of 
partons, similar to that in an ideal gas of gluons and quarks with Nf = 2. 

Using the same time evolution as given by Eqs. (|28|) and (0), we have evaluated 
the thermal dilepton yield with the above initial parton conditions, and the results are 
shown in Fig.^. Because of increased number of quarks and antiquarks, the contribution 
from the LO qq process is increased by a factor of about 16, while gluonic processes 
become less important in comparison. As a result, the total thermal dilepton yield 
is about a factor of six higher than that in the default scenario. Depending on the 
dilepton mass, this is still about a factor of 5 or more below the Drell-Yan yield. 

5.4 Dependence on initial parton conditions 

The thermal dilepton yield is sensitive to the initial conditions, such as tq and Tq, 
for the partonic matter. For example, the contribution from qq annihilation given in 
Eq.(|T7D depends on temperature as 

^ oc T^/^e-*^/^. (34) 

dM ^ ' 

The thermal dilepton yield can be large if the partonic matter starts at a very high 
initial temperature Tq. In the following, we study the dependence of the thermal 
dilepton yield on the parameters A^^O) \fi-, '^o and Tq. 

For a fully-equilibrated QGP at Tq = 570 MeV with iVy = 2 and an initial proper 
time To = 0.7 fm/c, the rapidity density of transverse energy is Q 

"^^^ roTTi^i^ f A + ^ll\ ^^T^ ^ 12, 000 GeV (35) 



dy " ^4 V15 20 

at central rapidity. This is more than an order of magnitude larger than most model 
estimates of the final dE^/dy. 

On the other hand, Eq.(|l|) with Xg^ = 0.060, Xq^ = 0.0072, tq = 0.7 fm/c and 
To = 570 MeV in the default scenario gives at central rapidity 

f?^ = f^?^.350GeV. (36) 

dy dy 4: 



Bose-Einstein and Fermi-Dirac distributions are used. 



Including additional Et from subsequent string fragmentation gives a final dE^/dy ~ 
900 GeV at central rapidity, which is comparable with predictions from most event 
generators and cascade models for heavy ion collisions at RHIC energies. 

In the following, we shall study the dependence of thermal dilepton yields on the 
initial conditions that have same value of initial dE'^/dy at central rapidity. According 
to Eg. (p9D , this is equivalent to taking initial conditions with the constraint 

I ^ diXi^o j tqTq = const. (37) 

We have already studied the effect due to the increase of the quark fraction in the 
partonic matter. Here, we study the dependence on initial conditions with a fixed 
Xgfi/Xgfi, i.e., with same values of both Ag^o^o^o^ and A^^o'^'o^o^- 

Using the same time evolution given by Eqs. (pSD and (pOf), the LO qq annihilation 
process gives 

dNgg ttR\ rt dRqq, , ri , ,or^.. fM\ 



iy-i5^^^-iy-^-'<^i^, 



dM^dY 2M Jro dM^ ' Jro " \T J 

ot ('^g,o''"o^o ) " dxx^Ki{x) (X / dxx^Ki{x). (38) 



3^0 3^0 



It is seen that this contribution depends only on one dimensionless variable M/Tq. 

In Fig.H, we show the value of the last integral in the above equation as a function 
of M/Tq. As Tq increases, the dilepton yield at a fixed mass M increases as the variable 
M/Tq shifts to the left, and the yield saturates at the high temperature limit. The slow 
increase and final saturation with increasing initial temperature Tq deviate from the 
naive expectation based on Eq. (|3^) . This is due to the constraint on the total dE^/dy 
in Eq.(^), which requires that as the initial temperature increases, the initial volume 
and/or the initial parton fugacities become smaller, thus making the total dilepton 
yield less sensitive to the initial conditions. We also see that the yield of higher-mass 
dileptons depend more strongly on the initial temperature. However, in mass range of 
2 — 4 GeV, where the yield of thermal dileptons is closer to the Drell-Yan yield, the 
variable M/Tq has a value between 3.5 to 7 when the default choice of Tq = 570 MeV 
is used. Fig.^ thus shows that the maximum increase one can obtain by increasing Tq 
while keeping the same total dE^/dy is only 15% at dilepton mass of 2 GeV and no 
more than a factor of 2 at dilepton mass of 4 GeV. Therefore, thermal dileptons are 
still well below the Drell-Yan yield even if we allow the maximum increase of the upper 
limit results in Fig. |^. 

Same studies can be carried out for the Compton and two-gluon processes, but they 
involve more parameters besides the dimensionless variable M/Tq due to the presence 
of the thermal quark mass fi. Relations similar to Eg . (pSf ) can still be obtained except 
for the logarithmic factors shown in Eq. (p2D . Similar conclusions thus hold as long as 
the screening mass is not extremely small. 
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6 Discussions and Summary 

In the present study, only leading-order processes have been considered. For low- 
mass dilepton production, higher-order effects have been studied |2^ using the Hard- 
Thermal-Loop (HTL) resummation method with effective propagators and vertices. It 
has been found that two-loop diagrams can give contributions at the same order as 
one-loop diagrams. Numerical calculations have shown that in heavy ion collisions at 
RHIC and LHC energies the two-loop contribution from the bremsstrahlung processes 
is important only for dileptons with mass M <^ T and becomes less important than 



the leading-order contribution for dileptons with mass M > 0.3 GeV [21|. Since our 
study mainly concerns dileptons with mass above 2 GeV {M ^ T), we do not expect 
that the inclusion of these processes will significantly change our results. 

However, there are still uncertainties in our results. First of all, the divergence in 
the dilepton production cross sections has been simply regularized by the thermal quark 
mass. Studies that take into account higher-order effects using the thermal field theory 
have indeed shown that for low mass dilepton or photon production in the QGP the 
thermal quark mass provides a self-consistent screening effect |jTl|, [l^, |T^ . However, the 
dependence of the dilepton production cross section on the thermal quark mass differs 
from that based on the lowest-order perturbative calculations as used in our study. 
Also, using different parton distributions can lead to different dependence of the cross 



sections on the thermal quark mass [|T^ . Since we have not carried out a self-consistent 
calculation including higher-order effects and have used the Boltzmann distribution 
instead of the quantum distribution, the coefficients that appear in association with 
the thermal quark mass as given in Eqs.(|ll]) and (|12D are expected to be different when 
a more complete calculation is carried out. To illustrate this uncertainty, we show in 
Fig.p the results from changing the screening mass by a factor of two but with the 
same default initial parton conditions. Since dilepton yields from the Compton and 
two-gluon processes roughly scale as log/i and log^/x, respectively, sizable differences 
are seen. 

We have only considered dilepton production from the stage where partons have 
reached thermal equilibrium. The contribution from the pre-thermal stage has not 
been included. Also, we have not considered chemical equilibration during the pre- 
thermal stage, which could change the initial conditions of the partonic matter such 
as the fugacities A^^o and Xg^. With the dilepton production cross sections given in 
Section |], the pre-thermal dilepton yield can be calculated using a pre-thermal parton 
distribution that is either parameterized or from a parton cascade code. With the 
latter, the entire time evolution of parton thermal and chemical equilibration can be 
studied, so the separation of pre-thermal and thermal stages, which is rather crude, 
is unnecessary. In this way, the evolution of parton fugacities can be self-consistently 
determined, and all secondary dileptons from the partonic stage can be calculated. In 
a previous study ||22[, we have carried out a calculation of secondary dileptons from 
the LO qq annihilation process within the parton cascade code ZPC ^3|, where initial 
parton fugacities were determined directly from the HIJING model [E3]. Since the 



default HIJING model gives too few quarks and antiquarks, with the ratio of the sum 
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of quarks and antiquarks to the total number of partons being only 4.5% (compared 
with the ratio of 15% in the default scenario of the present study based on minijet 
calculations using the perturbative QCD), the thermal dilepton yield in Ref.p2| is 
lower than the default result in this study by a factor of 10. 

In summary, we have studied thermal dilepton production in heavy ion collisions 
at RHIC energies. In particular, contributions from gluonic processes and their depen- 
dence on the initial parton conditions are investigated. We find that gluonic processes 
significantly enhance the thermal dilepton yield as one expects more gluons than quarks 
in the partonic matter formed in heavy ion collisions at RHIC. However, we have also 
shown that with the constraint of Eq.(^) on the total initial dE^/dy of the partonic 
matter, thermal dileptons are always significantly lower than those from the Drell-Yan 
process for any reasonable variations of the initial temperature, proper time, and fu- 
gacities. Our results show that thermal dileptons can compete with those from the 
Drell-Yan process only in a rather ideal case, when the initial energy density is higher 
than most current expectations, and parton chemical equilibration proceeds fast so 
that quarks are close to relative equilibrium with gluons. 

Acknowledgments: We thank B. Kampfer, W.L. van Neerven and S.M.H. Wong 
for helpful discussions. This work was supported in part by the National Science 
Foundation under Grant No. PHY-9870038, the Welch Foundation under Grant No. 
A-1358, and the Texas Advanced Research Project FY97-010366-068. 



12 



References 

[1] K. Kajantie, J. Kapusta, L. McLerran, A. Mekjian, Phys. Rev. D 34 (1986) 2746. 

[2] R. Baier, B. Pire, D. Schiff, Phys. Rev. D 38 (1988) 2814. 

[3] J. Cleymans, I. Dadic, Phys. Rev. D 47 (1993) 160. 

[4] E. Shuryak, Phys. Rev. Lett. 68 (1992) 3270. 

[5] T.S. Biro, E. van Doom, B. Mliller, M.H. Thoma, X.N. Wang, Phys. Rev. C 48 
(1993) 1275. 

[6] E. Shuryak and L. Xiong, Phys. Rev. Lett. 70 (1993) 2241. 

[7] B. Kampfer and O.P. Pavlenko, Z. Phys. C 62 (1994) 491. 

[8] M. Strickland, Phys. Lett. B 331 (1994) 245. 

[9] C.T. Traxler and M.H. Thoma, Phys. Rev. C 53 (1996) 1348. 

[10] X.-N. Wang, Phys. Rep. 280 (1997) 287. 

[11] T. Altherr and P.V. Ruuskanen, NucL Phys. B 380 (1992) 377. 

[12] M.H. Thoma and C.T. Traxler, Phys. Rev. D 56 (1997) 198. 

[13] R. Baier, M. Dirks, K. Redhch, D. Schiff, Phys. Rev. D 56 (1997) 2548. 

[14] P.J. Rijken, W.L. van Neerven, Phys. Rev. D 52 (1995) 149. 

[15] R. Hamberg, W.L. van Neerven, T. Matsuura, NucL Phys. B 359 (1991) 343; W.L. 
van Neerven, E.B. Zijlstra, Nucl. Phys. B 382 (1992) 11. 

[16] M. Gyulassy, Y. Pang, B. Zhang, NucL Phys. A 626 (1997) 999. 

[17] B. Zhang, CM. Ko, B.-A. Li, Z. Lin, preprint |nucl-th/9907"0T7| . 



[18] D.K. Srivastava, M.G. Mustafa and B. Miiller, Phys. Rev. C 56 (1997) 1064. 

[19] The difference between the values of our parton fugacities (Ag^o = 0.060, A^^o = 
0.0072) and those in ReLig (A° = 0.05, A° = 0.008) for the same parton number 
densities is due to the fact that we use the factorized Boltzmann distribution 
instead the factorized Bose-Einstein and Fermi-Dirac distributions. 

[20] P. Aurenche, F. Gehs, R. Kobes and H. Zaraket, Phys. Rev. D 58 (1998) 085003; 
iHd 60 (1999) 076002. 

[21] D. Pal and M.G. Mustafa, Phys. Rev. C 60 (1999) 034905. 



[22] S. Bass et al, preprint |nucl-th/9907090 

13 



[23] B. Zhang, Comp. Phys. Comm. 109 (1998) 193. 

[24] X.-N. Wang and M. Gyulassy, Phys. Rev. D 44 (1991) 3501; iMd 45 (1992) 844. 



14 



a) 



b) ^666^) 




N/x/x/v 




^^(^^^^ — ^ 



< K/\/\/v 




c) ^6^ — ^ 

'\/\/\/V 

^6^^ — 




'X/X/X/V 




-^6^ — >- 









Figure 1: Diagrams for dilepton production: a) qq annihilation, b) the Compton 
processes, and c) two-gluon fusion processes. Curly lines denote gluons, wavy lines 
denote virtual photons, and straight lines with arrows to the right (left) denote quarks 
(antiquarks). 
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Figure 2: Thermal dilepton rates from a quark-gluon plasma in full thermal and 
chemical equilibrium at a temperature T = 570 MeV. 
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Figure 3: Thermal dilepton rates from a gluon-dominated, under-saturated quark- 
gluon plasma with T = 570 MeV, \g = 0.060 and \g = 0.0072. 
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Figure 4: Thermal dilepton rates from a quark-gluon plasma in full thermal and 
chemical equilibrium at T = 570 MeV after being multiplied by the fugacity factors 



\g = 0.060 and A^ = 0.0072. 
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Figure 5: Time evolution of parton fugacities and temperature. 
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Figure 6: Thermal dilepton yields for central Au+Au collisions at RHIC with initial 
conditions of Tq = 570 MeV, Xg^ = 0.060, A^^o = 0.0072 and the time evolution 
according to Eqs. p8|) and (pOl). 
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Figure 7: Upper limit of thermal dilepton yields for central Au+Au collisions at RHIC 
with initial conditions of Tq = 570 MeV and A^^o = ^q,o = 0.028. 
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Figure 8: Dependence of the thermal dilepton yield from the LO qq annihilation on 
the initial temperature Tq with fixed rapidity density of the initial transverse energy. 
See Eq.(^) for the unit of the ordinate. 
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Figure 9: Screening mass dependence of thermal dilepton yields for central Au+Au 
collisions at RHIC with initial conditions of Tq = 570 MeV, A^.o = 0.060 and Ag^o = 
0.0072. Upper and lower curves correspond to results using half and twice the thermal 
quark mass as the screening mass, respectively. 
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